Pairing strength and symmetries of (K, Tl)Fej:Se2 in comparison with iron pnictides 
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We address singlet pairing states in the newly discovered superconducting iron selenides (K,Tl)FeiSe2 
within an incipient Mott approach. We analyze a five orbital t-Ji-J2 model, and contrast the results with those 
of a related model for the iron pnictides. Similar to the pnictides case, an Aig s^iy-i state and a Big d^2_y2 
state are quasi-degenerate in the Ji ~ J2 magnetic frustration regime. In contrast to the pnictides case, an 
Aig s^2^y2 state is not competitive in the relevant parameter regime. The dominant pairing channels are fully 
gapped. The pairing amplitudes are of the same order as their pnictides counterparts, making it natural that the 
maximum Tc is comparable in the two systems. 
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Introduction: An important characteristic of iron-based 
superconductors is their large materials phase space. Compar- 
ing the properties of the different materials families promise 
to shed considerable light on the physics of these systems. A 
major development along this direction is the very recent dis- 
covery of high temperature superconductivity in a new fam- 



I and (K, Tl)Fea;Sc2 
01 have a maximal 



ily of iron chalcogenides KFCj,Se2 
yfl. These and related compounds 
Tc ~ 30 A', which is comparable to that of 122 iron pnic- 
tides such as the doped BaFc2As2. Importantly, both the an- 
gle resolved photoemission experiments and ab initio 
band structure calculations using local-density approximation 
(LDA) li UTlll show that the fermiology of the new chalco- 
genides in the normal state is dramatically different from that 
in iron pnictides. In contrast to the pnictides where both elec- 
tron and hole pockets are present respectively at M and F 
points of the extended Brillouin zone, only electron pockets 
are present in the new chalcogenides. Equally important, the 
discovery has added a new twist to the consideration regarding 
the strength of the electronic correlations in the iron based su- 
perconductors. Similar to the iron pnictides, superconductiv- 
ity occurs near an antiferromagnetically ordered phase. How- 
ever in contrast to the pnictides where the antiferromagnetic 
phase is a semi-metal, the antiferroma gnet ic phase in the new 
iron chalcogenides is an insulator iB I12I1 . The observation 
of an insulating state and the emergence of antiferromagnetic 
order with large moment f\3\ in the absence of quasi-nested 
electron and hole pockets ceitainly point towards the impor- 
tance of correlation effects. It may be that the parent com- 
pounds of the new materials are Mott insulators, as suggested 
by microscopic model calculations lfl4i [Tsll that incorporate 
the band-narrowing effect seen in LDA studies lflo[[l6ll . 

Ever since the discovery of the high temperature supercon- 
ductivity in iron pnictides, the strength of the electronic cor- 
relation has been the subject of debate. Both weak and strong 
coupling pictures had partial success in explaining the (tt, 0) 
magnetic order in the parent compounds, and the emergence 
of a sign changing Aig state as the dominant pairing channel. 
In the parent iron pnictides, a large electrical resistivity ("bad 
metal"), a strong suppression of Drude weight lllTL as well as 
the temperature-induced spectral weight transfer ni8[ - l20ll sug- 
gest significant coiTelation effects. These properties have been 



interpreted in terms of a proximity to a putative Mott transition 
I 2ll - l23[l . Interestingly the room temperature electrical resis- 
tivity of the new iron chalcogenides is several orders of magni- 
tude larger than that of pnictides, and with decreasing temper- 
ature the resistivity shows exponential enhancement H, 12]. 
Strong correlation effects are also observedjii the suppression 
of the low frequency optical conductivity [i24<1 . All these sup- 
port the notion that both the parent iron selenides and parent 
iron pnictides are in the vicinity (but on two different sides) 
of the Mott transition. They make it important to understand 
the superconducting pairing in these materials from a strong 
coupling perspective 12114231 l25l - l29[l . 
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FIG. 1. (Color online) The contrasting Fermi surfaces of (a) 
KFea;Se2, and (b) iron pnictides, in the extended Brillouin zone cor- 
responding to one iron per unit cell, as obtained by using a five orbital 
tight binding model, for electron doping 5 — 0.15. There are only 
electron pockets at M points (±7r, 0) and (0, ±7r) for KFea;Se2. For 
iron pnictides, there are in addition two hole pockets at the F point 
(0,0). 

Several experimental results have especially motivated 
our study. The comparable maximum Tc observed in the 
(K, Tl)Fe3.Se2 compounds and their 122 pnictides counter- 
part suggest some degree of commonality in the underlying 
mechanism for superconductivity. At the same time, even 
though the pairing symmetry remains to be determined, nearly 
all experiments have pointed towards a fully gapped super- 
conductivity. This is in strong contrast to the iron pnictides, 
in which both nodeless and nodal features have been seen in 
different materials families or even at different dopings of a 
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same family. Our results will shed light on such similari- 
ties and differences. We note that weak coupling approaches, 
based on functional renormalization group BOiblll and spin- 
fluctuation exchange RPA l32l 13311 . have led to a dominant 



instability in the Big dj.2„j,2 pairing channel due to the repul- 
sive inter-electron-pocket pair scattering. Our work is closer 
in spirit to the strong-coupling approach to the (K, Tl)Fc2;Se2 
systems lfl4l[l5ll . although the pairing strength and aforemen- 
tioned similarities and differences compared to the iron pnic- 
tides have not been studied. 

Five-orbital Model: Our multi-orbital t-Ji-J2 Hamilto- 
nian |21, 22, 25l|23l is given by 



H 



i<j,a..(3,s i,a 



^ J"^ (Sia ■ SjfS - -UialljA 



(1) 



where c\^^ creates an electron at site i. 



and spin projection 
the hopping matrix. 



with orbital index a 



3; ^ is the chemical potential and t' 
The orbital index a = 1,2, 3, 4, 5 re- 



m 



spectively correspond to five 3d orbitals 3d2;^, 3dyz, 3d, 
"iAxy, and 3d3;j2_^2 of iron. The nearest-neighbor (n.n. 
and next-nearest-neighbor (n.n.n., {{ij))) exchange interac- 
tions are respectively denoted by Ji 
operator S^a = 5 E...' 4c.s^s, 



and . The spin 
and the density oper- 
ator riia = J2s Cias'^ias, with (7 representing the Pauli ma- 
trices. The double-occupancy prohibiting constraint from the 
fermion is implicitly incorporated by the renormalization of 
the dispersion. As discussed at length in the Supplementary 
Material ll34ll . the phase diagram is insensitive to this when the 
exchange interactions are measured w.rt. the renormalized 
band dispersion. In a previous study of a multi-orbital t-Ji- 
J2 model for iron pnictides 1271], it has been demonstrated by 
three of the present authors that the dominant pairing symme- 
try is governed by the intra-orbital exchange interactions, and 
the consideration of the orbitally off-diagonal exchange inter- 
action only introduces quantitative modifications of the phase 
diagram. Therefore to keep the analysis simple, we will con- 



sider the orbitally diagonal exchange couplings J"'^ ^ JiSa 



and = J2'5q^. Experiments have suggested that the or- 
dering of the Fe-site vacancies goes away in some supercon- 
ducting materials 35, 3^, but not in others H, 13|. For 
simplicity, we will consider here the case with the 122 tetrag- 
onal symmetry. 

To describe the fermiology, we use a tetragonal symmetry 
preserving tight binding model, involving all five 3d orbitals 
of iron. In the momentum space the 5x5 tight-binding ma- 
trix, and its eigenvalues will be respectively denoted by £^k, 
and Ej^]f^. The total number of electrons is determined by the 
chemical potential and the dispersion relations, according to 
n ~ 2 J2j k ^(m ^ -^j,k)- and carrier doping 5 = |n — 6|. The 
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FIG. 2. (Color online) Zero temperature phase diagrams of (a) 
KFea;Se2 and (b) iron pnictides, for electron doping S = 0.15. Re- 
gion I corresponds to Aig state with dominant 5^.2^2 channel. Re- 
gions II and III mark an Aig + iBig state with dominant s^2y2 



and d_. 



channels (II) and dominant s 



and d_ 



chan- 



nels (III); the phase locking occurs only at low temperatures. Region 
IV for KFea;Se2 is a pure Big 



state with dominant d^2_y2 channel. 
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FIG. 3. (Color online) The pairing amplitudes for xzlyz orbitals 
of KFeiSe2 for electron doping 5 = 0.15. In the J2 dominated 
regimes, for J2 /D < 0.3, Aig s^2y2 is the strongest pairing channel. 
In the likewise Ji dominated regime, this is the Big d^2_y2 channel. 



tight-binding parameters are fixed by fitting the band struc- 
ture obtained from the LDA calculation. The details of the 
tight-binding parametrization for iron chalcogenides and the 
associated band dispersions are discussed in the Supplemen- 
tary Material ||34|| . For (K, Tl)Fca;Sc2 the band dispersions 
correctly produce the electron pockets near the M points, as 
illustrated in Fig. |l(a)| for KFcj:Se2. The electron and hole 
like Fermi pockets obtained from a similar five orbital tight- 
binding model of iron pnictides are shown in Fig. |l(b)| to 
contrast the fermiology of the two materials. In the Supple- 
mentary Material ||34|| we show the fermiology of TlFea;Se2, 
which again consists of only electron pockets. We note that 
photoemission experiments jsj-lztl have suggested that very 
weak electron-like pockets also exist near the F point; unlike 
their hole counterpart in the pnictides, these electron pockets 
have very small spectral weight and are not expected to play 
an important role in driving superconductivity. 
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Degenerate spin singlet pairing states in the absence of ki- 
netic energy: The spin singlet, intraorbital pairing operators 
are defined as A^^^a = {cra^c,+eai - Ci„iCi+eat)/2, where 
e = x,y,x ± y. Two types of pairing states q, = ±Ay 
respectively denoted as Sx2_^_y2 and d^2_y2 states with mo- 
mentum space form factors gx^+y^.k = coskx + cos ky, 
gx^-y^^k = cosfca; — cos ky arise due to the nearest neigh- 
bor exchange interaction, and they are energetically degen- 
erate in the absence of kinetic energy. Similarly the next 
nearest neighbor exchange gives rise to two degenerate pair- 
ing states Ax+y,a = ^Ax~y.a, respectively denoted as Sx2y2 
and dxy states, with momentum space form factors gx^y^^k = 
cos{kx — ky) + cos{kx + ky),gxy^k = cos(kx — ky) + cos{kx + 
ky). A strong magnetic frustration, characterized by Ji ^ J2, 
leads to an enhanced degeneracy among the four paired states. 

Degeneracy lifting effects of the kinetic energy: The 
kinetic energy term lifts most of these degeneracies and 
in the strong frustration regime leaves a quasi-degeneracy 
among a reduced set of pairing states. To study the 
full problem we perform a mean-field decoupling lEsl 



0.15 



27L 13711 of the exchange interaction terms. We intro- 
duce four complex singlet pairing amplitudes for each 
orbitals, and write the following 5x5 pairing matrix 
Ak = Eadiag[A« n,A« 22,A« 33,A« 44,A« 55], where 
aa ^ ^QQffa.ki ^nd the indcx a corresponds to 
Sx2j^y2, dx2_y2, Sx2y2 and dxy symmetries. In the sub- 
space of xz and yz 3d-electron-orbitals, which transform 
as a doublet under the tetragonal point group operations, 
there are the following four classes of intra-orbital pair- 
ing states for an orbitally diagonal J1-J2 model; (i) Aig : 

[Sx2ly2gx^+y^M + st2y2gx^y^,k\T0 + d^2'_y2gx-^ -y^ M^z; &) 



Bi 



2gx^ 



Alg 1 

Sx2y2gx^y^,k\T0 

r -Big 

l^x'^+y^9x^+y^M- 



2.k7"0- 



^x^y^y^^y 



dxy" gxy.kTf). Here tq 



(iii) A2g : dxy" gxy.kTz\ and (iv) B2g 
and Tz are the Pauli matrices in the two-orbital subspace. The 
eight pairing amplitudes s^2^j^2 etc. are obtained as linear 
combinations of AJj^ and A22, which are in turn linear com- 
binations of Ae.ll, and Ae.22- 

Pairing phase diagram: We determine the pairing gap 
matrix by minimizing the ground state energy density 



/ 



£^k,j+M), (2) 



with respect to all Ae,tM, where ±fk,i are the eigenvalues of 
the mean-field Hamiltonian 



Ck - A«l5x5 



Ak 

~Ck + Ml5x5 



*k, (3) 



and 'I'jj. = (cjj.^,|., C-kai) is the ten-component Nambu spinor. 
Our calculations are performed in a canonical ensemble for a 
fixed doping or electron density. The grand canonical results 
for the phase diagram are qualitatively similar to the canonical 
ensemble results. 

We show two illustrative zero temperature phase diagrams 
in Fig. |2(a)| and Fig. |2(b)| respectively for KFca;Sc2 and iron 
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FIG. 4. (Color online) Comparison of the competing dominant pair- 
ing amplitudes Aig 3.^,2^2, Aig 3^2+^2, Big d^2_y2 for xzlyz or- 
bitals of (a) KFe^Se2 and (b) iron pnictides, both for electron doping 
5 = 0.15. For KFeiSe2, the amplitude for the Aig s^2^y2 channel 
is strongly suppressed compared to the pnictides case. Correspond- 



ingly, a pure Bij 
in iron pnictides. 



2 state is observed for KFea;Se2 but is absent 



pnictides, with an electron doping 5 = 0.15 (n = 6.15). The 
phase diagrams are shown for < Ji, J2 < 0.3-D, where D 
is the band width. The phase diagram in Fig. |2(a)| is separated 
into four regions I, II, III, and IV, whereas that in Fig. |2(b)| is 
separated into three regions I, II, and III. In the J2 dominated 
region I, a pure Aig state occurs with Sx2y2 being the dom- 
inant component. For iron pnictides this state corresponds 
to the s± state, that changes sign between electron and hole 
pockets. Close to the Ji ^ J2 magnetic frustration regime, 
the Aig Sx2y2 and Big dx2^y2 states are quasi-degenerate. 
In region II, due to this quasi-degeneracy between two types 
of pairing symmetries which follow different crystallographic 
representations, the lowest energy paired state is a time rever- 
sal symmetry breaking Aig + iBig state. In the Ji dominated 
regime III, the paired state is again a time reversal-symmetry 
breaking Aig + iBig state, but the dominant components are 
Aig Sj.2_,^j,2, and Big dx2^y2 pairings. Note that there is only 
a crossover between regions II and III, and time-reversal sym- 
metry breaking is expected only at very low temperatures. In 
the phase diagram of Fig. |2(a)| for KFe:i.Se2, the additional 
region IV is dominated by Ji, but in this region a pure Big 
d2:2_j^2 state occurs. Notice that this pure Big phase is absent 
for iron pnictides. To illustrate the robustness and applicabil- 
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ity of the phase diagram Fig. |2(a)| for KFe2;Se2, we show a 
quaUtatively similar phase diagram for TlFc3;Sc2 in the Sup- 
plementary Material [34]. 

To get an overall picture of the pairing strength, we show 
in Fig. [3] a three dimensional plot of the pairing amplitudes 
for xzlyz 3d-electron-orbitals of KFe2:Se2. The other three d- 
orbitals may also contain significant pairing strength, depend- 
ing on their spectral weight on the Fermi pockets (see ||34|] for 
details). Similar results for iron pnictides are also shown in 
the Supplementary Material ll34ll . The competition among the 
dominant pairing channels can be clearly observed in Fig. |4(a)| 
and Fig. |4(b)| respectively for (K, Tl)Fe2:Se2 and iron pnic- 
tides. In these figures we show the dominant pairing ampli- 
tudes as a function of Ji/D, for a fixed value of J2/D = 0.1. 
In Fig. |4(b)[ we notice, apart from the absence of pure Big 
phase, that Aig Sj;2^y2 state has a significant contribution in 
the entire regime < Ji < OAD. This is important to under- 
standing the nodeless vs. nodal gaps in the Aig state of iron 
pnictides. This pairing state by itself will have nodes on the 
electron pocket, and for KFc2;Sc2 it is never very competitive 
in the weak to moderate Ji regime. Corresponding pairing 
amplitudes of TlFea;Sc2 are shown in the Supplementary Ma- 
terial IMI- 

The role of the fermiology behind the emergence of the 
dominant pairing states can be understood by considering the 
strength of the pairing kernels for different channels. Due to 
the smallness of dxy form factor sin kx sin ky, in the vicinity 
of M and F points, the dxy pairing kernel is very small for 
a set of small Fermi pockets, relevant to the iron based su- 
perconductors. For (K, Tl)Fe2;Se2, which only possess small 
electron pockets at M points, the pairing kernel for Sx2+y2 is 
also very small. In contrast both d^2 _y2 and Sx2y2 pairing ker- 
nels are much bigger around the M points, and in the weak to 
moderate ( Ji, J2) regime, these two states emerge as the dom- 
inant pairing states. For iron pnictides there are two small 
hole pockets around F point, and due to these pockets Sx2y2 
kernel becomes larger than dx2-y2, and s^2j,2 dominates the 
larger portion of the phase diagram. Again due to the pres- 
ence of the hole pockets Sx2^y2 has a stronger competition 
with dj,2_,y2 state in pnictides. 

Experimental implications: For (K, Tl)FexSc2, our cal- 
culations lead to dominating Aig Sx^y^ and Big dx2-y2 states 
in competition. Since both have nearly isotropic and node- 
less gaps on the electron pockets, our results are consistent 
with existing experiments. For iron pnictides, the dominance 
of Aig Sx2y2 state in a large portion of the phase diagram 
is consistent with the expected s± paired state predicted by 
both weak and strong coupling approaches. The sizable con- 
tribution from Sx2^y2 channel to Aig state is important to un- 
derstanding a possible strong momentum dependence or gap 
nodes of the Aig state on the electron pockets in different 
pnictides compounds. 

At the same time, our calculation shows similar pairing am- 
plitudes for the two systems, thereby giving a natural under- 
standing for the comparable transition temperatures. 

At present there is no experimental estimation of the ex- 



change couplings in the paramagnetic normal state. Future 
inelastic neutron scattering experiments should provide guid- 
ance in such a direction, and this will be important for under- 
standing the pairing symmetry of (K, Tl)Fc2.Sc2 compounds. 
We close by noting that recent experiments have suggested 
coexistence of a large moment antiferromagnetic order and 
superconductivity |i3| or a phase separation Theoreti- 
cally, how superconducting and magnetic orders interplay in 
these materials is an intriguing question for future studies. 
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SUPPLEMENTARY MATERIAL FOR EPAPS 



I. DETAILS OF TIGHT-BINDING PARAMETRIZATION 



For the five orbital, tetragonal symmetry preserving tight binding model, we adopt the parametrization of Refl38l We have 



fitted the LDA band structure with the band dispersion found from the tight-binding model, to determine the hopping parameters. 
For iron pnictides we use the tight binding parameters of Ref|38l For (K, Tl)Fe2;Sc2, we have performed calculations with two 
different sets of hopping parameters. These two sets respectively were derived from fitting the LDA results for KFe2Sc2, and 
TlFc2Se2. In the main text we have shown the electron pockets and pairing phase diagram for the band structure corresponding 
to KFcx Sc2 . In this section we provide the explicit parametrization and band dispersions for both tight-binding models, and also 
show the electron pockets derived from TlFc2Se2. 




FIG. 5. The band dispersions of (a) KFe2Se2 and (b) TlFe2Se2, along the high-symmetry directions of the extended Brillouin zone (one iron/ 
unit cell). The band dispersions have been obtained by diagonalizing the appropriate five orbital tight binding model. In panel (c) we show the 
Fermi surfaces of TlFeiSe2 with electron doping S — 0.15, which consist of only electron pockets at M points. 
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Tight-binding model for KFe2Se2." The following 5x5 tight binding matrix for KFea;Se2 has been used. The orbital 
indices 1, 2, 3, 4 and 5 respectively correspond to dxz, dyz, d^2_y2, d^y, and d^z^-r'^ orbitals. 

= -0.36559- 0.1095 cos fc^ - 0.81736 cos fey + 0.83524 cos /c^r cosky + 0.03114(cos2A::r - cos2fcy) 

-0.03464 cos 2fc^ cos fey - 0.12572 cos 2fey cosfc^ + 0.07596 cos 2fc^ cos2fey 
= -0.36559 - 0.81736 cos fe^ - 0.1095 cos fey + 0.83524 cos fe^ cosfey - 0.03114(cos2fe^ -cos2fcy) 

-0.12572 cos2kxCosky - 0.03464 cos 2fey cos k^ + 0.07596 cos 2fca; cos 2fey 
= -0.56466 + 0.65046(cos k^ + cos fey) - 0.39132 cos k^ cos fey - 0.01074(cos 2kx + cos 2ky) 
Ck** = -0.00096 + 0.41266(cos k^ + cos fey) + 0.26328 cos k^ cos fey - 0.0705(cos 2fe^ + cos 2ky) 

— 0.08756(cos2fca; cos ky + cos 2fey cos k^) + 0.01692 cos2fca; cos 2ky 
£,1^ = -0.91583- 0.0854(cosfc^ + cos fey) + 0.02234(cos2fc^ + cos2fey) 

+0.00708(cos2fca; cosky + cos2fey cosfe^;) — 0.05396 cos 2fca; cos2fcy 
^k^ = -0.40644sinfc^sinfey + 0.08068(sin2fc^sinfcy + sin 2 fey sin fe^;) - 0.13092 sin2fe:j; sin2fey = ^k^* 
^k^ = -0.62894zsinfey + 0.249isinfeycosfe;j; - 0.0412z(sin2fey cos fc^; - cos 2 fc^: sin fey) = £,1^* 
^23 ^ 0.62894i sinfe^, - 0.249z sin fc^; cosfey + 0.0412z(sin2fc^ cosfey - cos 2fey sin fe^;) = £,1^* 
Ck** = 0.2757isinfe;j; + 0.0042i cosfcj^ sinfc^^ + 0.0416i sin2fc:j, cosfey = £l^* 
£1^^ = 0.2757isinfey + 0.0042icosfc:j,sinfcy + 0.0416i sin2fcy cos fe-^, = £^^* 
£l^ = -0.0965isinfcy + 0.40384isinfeycosfca; +«(0.04816)sin2fcycos2fe:r = £l^* 
^25 ^ -0.0965isinfc^ + 0.40384isinfc^cosfcy + 0.04816i sin2fe:r cos2fey = £l^* 
^ _0.1918(sin2fey sinfc^r - sin2fc^sinfcy) = £l^ 

= -0.61932(cosfca; - cosfey) - 0.05992(cos2fc3; cosfey - cos2fey cosfey;) = ^k^ 
£t^ = -0.33436 sin fc^ sin fey - 0.03064 sin 2fe^ sin2fcy = £l'^ (4) 



TlFe2Se2 .• The following 5x5 tight binding matrix £k has been deriving from fitting the LDA results of TlFc2Se2. 

^k^ = -0.35956 - 0.48396 cos k^ - 0.69426 cos ky + 0.93156 cos k.^ cos ky + 0.20428(cos 2k^ - cos 2fey ) 

-0.23556 cos 2fe3: cosfey + 0.2714 cos 2fey cosfe-^; + 0.02912 cos 2fe-:r cos2fey 
Ck^ = -0.35956 - 0.69426 cos k^ - 0.48396 cos ky + 0.93156 cos fc^ cos ky - 0.20428(cos 2fc^ - cos 2fey) 

+0.2714 cos 2fca; cos fey - 0.23556 cos 2fey cosfc^^ + 0.02912 cos 2fc:r cos2fey 
£l^ = -1.11574 + 0.76802(cosfe;j; + cos fey) -0.10188 cosfex cos fey + 0.02612(cos2fc^ + cos 2fey) 
£t'^ = 0.09324+ 0.76338(cosfe^ + cosfej^) + 0.63348 cosfey cosfey - 0.06164(cos2fc^ + cos2fey) 

— 0.18652 (cos 2fca; cos fey + cos 2fey cos fc^;) — 0.1074 cos2fc2: cos 2fcy 
£l^ = -0.74545 + 0.07376(cos fe^ + cos fey) -0.05124(cos2fe^ + cos 2fcy) 

— 0.12872(cos2fc^ cos fey + cos 2fey cos fc^) + 0.08408 cos2fc^ cos 2fey 
£l^ = -0.97632 sinfc^ sin fey - 0.01744(sin2fej; sin fey + sin 2 fey sin fc^) - 0.04692 sin 2fe2; sin 2fey = ^k^* 
^k^ = -0.73408isinfey - 0.08848i sin fey cosfey; + 0.10508i(sin2fey cos fe^; - cos2fca; sinfcy) = £l^* 
£l^ = 0.73408isinfc2: + 0.08848z sinfe^; cosfey - 0.10508i(sin2fca; cos fey - cos2fey sinfe^^) = ^k^* 
Ck* = 0.46354isinfc:j, - 0.22408i cosfey sin fe^, + 0.15492i sin2fe^ cos fey = £^^* 
£l^ = 0.46354isinfcy - 0.22408i cosfcx sin fe-y + 0.15492i sin2fcy cos fe;^; = Ck^* 
£l^ = -0.26344isinfey + 0.24328i sinfcy cos fe;^; - 0.11724i sin2fcy cos 2fc:j, = £l^* 
£l^ = -0.26344isinfe3; + 0.24328i sinfe^; cosfey - 0.11724* sin 2 fc^ cos 2 fey = ^k^* 
^k" = -0.20372(sin2fey sinfe^; - sin2fe2; sin fey) = ^k^ 

£l^ = -0.4617(cosfca; - cosfey) + 0.13144(cos2fca; cosfey - cos2feyCosfc^) = £l^ 

^45 ^ -0.6828 sin fc^ sin fey - 0.1088 sin 2fc^ sin 2fcy = £l'^ (5) 
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II. PAIRING AMPLITUDES FOR TlFexSea AND IRON PNICTIDES 

In this section we first show the pairing phase diagram TlFe2^Se2 in Fig. |6(a)[ for electron doping 5 ~ 0.15. Also shown, 
Fig. |6(b)| are the dominant pairing amplitudes projected onto the xzlyz 3d-orbitals. The phase diagram and the strength of the 
pairing amplitudes are similar to those of KFc2;Sc2, described in the main text. In Fig. |7(a)| and Fig. |7(b)[ we respectively show 
three dimensional plots of the pairing amplitudes for TlFea;Se2 and iron pnictides, both with an electron doping 5 = 0.15. 




FIG. 6. In panel (a) we show the zero temperature phase diagram of TlFea;Se2 for an electron doping 5 — 0.15. The phase diagram is 
similar to that of KFea;Se2 (Fig. |2(a)] (. The regions I, II, III, and IV respectively correspond to an Aig state with 5^.2^2 as the dominant pairing 
channel, a time reversal symmetry breaking Aij + iBig state with s^2y2 and d^2_y2 as the dominant Aig and Big pairing channels, a likewise 
Aig + iBig state with s^2_^_y2 and d^2_y2 as the dominant Aig and Big pairing channels, and a pure Big state with d^2_y2 pairing channel. 
In panel (b) we show the competing dominant pairing amplitudes Aig s^2y2 , Aig s^2_^_y2 , Big s^2y2 for xzldy^ orbitals of TlFea;Se2 with an 
electron doping 5 = 0.15, and J2/D = 0.1. 




(a) (b) 



FIG. 7. The pairing amplitudes for xzlyz orbitals of (a) TlFej:Se2 and (b) iron pnictides with electron doping 5 = 0.15. Compared to 
TlFei,Se2, iron pnictides have stronger Aig s^2^y2 pairing. 



III. PAIRING PHASE DIAGRAMS OF A TWO ORBITAL MODEL WITH AND WITHOUT DOPING DEPENDENT BAND 

RENORMALIZATION 

To demonstrate the robustness of our pairing phase diagram against the doping dependent band renormalization effects, we 
consider a two orbital model involving only xz and yz orbitals. For simplicity we consider hole doping; the electron doping case 
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can be treated in a similar manner after performing a particle-hole transformation. The Hamiltonian of interest is 



i<j ,<y,^ ,a i,a.,(T i'Cj,a.f3 



(6) 



with the double occupancy prohibiting constraint J2a '^aa'^iaa- ^ 1 for each orbital. The constraint is imposed by introducing 
slave boson operator bi and fermionic spinon operator fi for each orbital ll37ll and the Hamiltonian is transformed to 



H 



i<j.aJ3 .(T i,a.,{T 



E ■ 

■i<j,a,l3 



2 



(7) 



where Bj^ = {fja-\-.fjpi ~ /Li/J/at) '■^^ spin-singlet pairing operator for the fermionic spinous, and Ai.a's are Lagnrage 
multipliers which enforce the occupancy constraints. At zero temperature the slave bosons are bose condensed and the bose 
operators have finite expectation values. In the tetragonal symmetry breaking particle-hole channel order, we have = 
(^1,2) = \/5/2, where S is the hole doping and {Xi^i) = {Xi,2) = A, which is absorbed into the chemical potential. The 
expectation value of the bose operators renormalizes the kinetic energy term by a factor of 6/2. As in the main text we choose 
the interactions to be diagonal in the orbital space. After mean-field decoupling, the free energy density 



6 S 
k+ + fk- - ^E^+ - ^Eu- + 2/iX8) 



is minimized with respect to pairing amplitudes with the constraint 711 = 712 = 1 — |. The quasiparticle dispersions fki in the 
paired state are calculated from the 4x4 Nambu matrix 

where 



56 
2 



Ak 

■ + A*l2x2 



(9) 



Ji(A cos(fc^) + A 



^x — y.aa COs(^kx ^y)) 



(10) 



For the band structure we choose the two orbital tight-binding model of Ref. |39l The 2x2 tight-binding matrix ^k = = 
^k+l2x2 + S.'k-Tz + ^UxyTx, where Pauli matrices r^) operate on the orbital indices, and ^k+ = —{ti+ i2)(cos + cos ky) — 
4^3 cos kr^ cos ky, ^k- = ^(^1 ^ t2)(cos k^ — cos ky), ^ka;y = —4^4 sin kj. sin fey are respectively Aig, Big and 62^ functions. 



The band dispersion relations i?k± = Ck+ ± y ^k- + ^kxy^ 8^^^ electron pockets at k = (tt, 0) and (0, tt), and two 

hole pockets at k = (0, 0) and (tt, tt). The following values of the hopping parameters, ti = —t, t2 = l.it, ^3 = ^4 = — 0.85t 
were obtained in Ref.l39l by a fitting of the LDA bands. The quasiparticle dispersions in the paired state are given by 



■k± 



''5Ck+ 



p 2 2 

(Ck- + Ckxa 



lA 



k,ll| 



lA 



k,22| 



lA 



k,ll| 



lA 



k,22| 



\z.xy 



^|Ak,ll — Ak,22 




(11) 



After explicitly accounting for the occupancy constraints and associated band renormalization we obtain the pairing phase 
diagram shown in Fig. [8(a) A qualitatively similar phase diagram for this model, obtained by ignoring the band renormalization 
effects, is shown in Fig. 8(b)| 



IV. ORBITAL CHARACTER OF THE FERMI SURFACE AND PAIRING AMPLITUDES FOR xy ORBITAL 



In this section we first show the orbital weights on the electron pockets for both KFe^Se2 and pnictides, for electron doping 
6 = 0.15 respectively in Fig. |9(a)| and Fig. |9(b)| For pnictides, the xzlyz orbitals have the dominant orbital weights on the 
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FIG. 8. Qualitatively similar zero temperature phase diagrams of a two orbital t — Ji — J2 model obtained (a) by considering a doping dependent 
renormalization of the band dispersions due to the local constraint of zero double occupancies and (b) by ignoring the doping dependent band 
renormalization. In (a) the effective hopping parameters U — ti^. The regions I, II, and III respectively correspond to an Aig state with s^2y2 
as the dominant pairing channel, a time reversal symmetry breaking Ai^ + iBig state with s^2y2 and d^2_y2 as the dominant Aig and Big 
pairing channels, and a likewise Aig + iBij state with 5^24.^2 and d-^2_y2 as the dominant Aig and Bij pairing channels. 




FIG. 9. The orbital weights on the electron pocket centered at (tTjO) of (a) KFea:Se2 and (b) iron pnictides for electron doping 5 = 0.15. 
6 is the winding angle of the pocket with respect to its center. The weights on the electron pocket centered at (0, tt) can be obtained by 
interchanging the xz and yz components and shifting 9hy n /2. 

two hole pockets (not shown), and also contribute significantly on the electron pockets. Whereas for KFc2;Se2 (and also for 
TlFG2;Se2), there are no hole pockets, and the contribution from xy orbital is considerably enhanced on the electron pockets. 
This indicates that the xy orbital plays a more important role in building the electron pocket of (K, Tl)Fe:cSc2 than in pnictides. 

The enhanced xy orbital character on Fermi surface of KFea;Se2 affects the electron pairing. To see this we show the pairing 
amplitudes at J2/D = 0.1 of a;y orbital for both KFc2^Se2 and pnictides at electron doping 5 = 0.15 in Fig. |10(a)| and Fig. [l0(b)| 
respectively. Comparing with Fig. |4(a)| and Fig. |4(b)| we see that for pnictides, the contribution to the three competing dominant 
pairings from the xzlyz orbital is comparable with the one from xy orbital. But for KFea;Sc2, the contribution from xy orbital 
can be stronger This is especially true in the regime Ji/ D < 0.07, where only Aig pairing is present. 

Overall, however, the pairing amplitudes for the xy orbital are similar between (K, Tl)FexSc2 and iron pnictides. This is 
seen in the three-dimensional plots given in Figs. |l l(a)| and [TT(b)| which are similar to their counterparts for the xzlyz orbitals 
(Figs.|3]and[7(b)). 
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FIG. 11. The pairing amplitudes for the xy orbital of (a) KFea;Se2 and (b) iron pnictides, both for electron doping 5 = 0.15. 



